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Under sufficient permanent random covalent bonding, a fluid of atoms or small molecules is 
transformed into an amorphous solid network. Being amorphous, local structural properties in such 
networks vary across the sample. A natural order parameter, resulting from a statistical-mechanical 
approach, captures information concerning this heterogeneity via a certain joint probability distribu- 
tion. This joint probability distribution describes the variations in the positional and orientational 
localization of the particles, reflecting the random environments experienced by them, as well as 
further information characterizing the thermal motion of particles. A complete solution, valid in 
the vicinity of the amorphous solidification transition, is constructed essentially analytically for 
the amorphous solid order parameter, in the context of the random network model and approach 
introduced by Goldbart and Zippelius [Europhys. Lett. 27, 599 (1994)]. Knowledge of this order 
parameter allows us to draw certain conclusions about the stucture and heterogeneity of randomly 
^ ■ covalently bonded atomic or molecular network solids in the vicinity of the amorphous solidification 

transition. Inter alia, the positional aspects of particle localization are established to have precisely 
the structure obtained perviously in the context of vulcanized media, and results are found for the 
analogue of the spin glass order parameter describing the orientational freezing of the bonds between 
particles. 
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I. INTRODUCTION AND OVERVIEW 



The purpose of this Paper is to address the statistical structure of the amorphous solid state via a simple model of 
■ a three-dimensional vitreous medium consisting of covalently bonded atoms (or low-molecular- weight molecules) [Q. 
We shall do so essentially analytically by making use of techniques drawn from the field of the statistical mechanics of 
systems having quenched randomness J2j . The model of vitreous media on which we shall focus is that introduced by 
Goldbart and Zippelius [gj, which takes as ingredients a thermodynamically large number of particles between which 
some large number of permanent random covalent bonds are introduced. The quenched randomness is encoded in the 
information describing which pairs of particles are covalently bonded; the remaining (annealed) degrees of freedom 
correspond to the unconstrained positions of the particles and the orientations of the orbitals. This model exhibits 
a continuous equilibrium phase transition from the liquid state to the amorphous solid state when the density of 
C3 introduced bonds exceeds a certain critical value. It is on the structure and heterogeneity of this state that we hope 
I ■ to shed some light. 

As an example of the type of medium we have in mind, consider networks formed by the polycondensation of Si(OH)4 
molecules, during which H2O is eliminated between pairs of hydroxyl (OH) groups on certain randomly selected pairs 
of Si(OH)4 molecules so as to form Si-O-Si bonds. The amorphous solidification of such media has been studied in 
many experiments; we cite as an example those of Gauthier-Manuel et al. [Q. As it is our intention to develop a rather 
general model of random networks, and to focus on universal properties, it is not necessary for us to incorporate the 
specific details of the medium. For example, we shall not be accounting for the bond geometry associated with the 
so-called bridging oxygen atoms between the silicon atoms. In the model, both types of Si orbitals, those connected 
to hydroxyl groups and those connected to bridging oxygen atoms, will simply be referred to as "orbitals'.' A second 
example of the type of media we have in mind is provided by amorphous silicon networks |IJ] , especially those in which 
hydrogen passivates bonds unconnected to other silicon atoms. 

The structural characterization of the vitreous state that we shall construct will be statistical in nature, reflecting 
the intrinsic heterogeneity of the environments that the constituent particles in vitreous media inhabit. It will take 
the form of a joint probability distribution characterizing the fraction of particles that are localized in the vitreous 
state, and will describe the spatial extent of the thermal fluctuations in their positions, the degree and character of 
the thermal fluctuations in the orientations of the orbitals that are capable of participating in covalent bonds, and 
the strength and nature of the correlations between the thermal fluctuations in the particle positions and the orbital 
orientations. Moreover, rather than dealing with media having a specific architecture (i.e., a specific realization of 
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introduced bonds), we shall consider an ensemble of architectures, all characterized by a common parameter governing 
the probability that a permanent chemical bond was formed between any pair of nearby orbitals. 

A statistical description of an amorphous solid state has previously been developed and explored in the context of 
vulcanized (i.e., randomly permanently crosslinked) macromolecular media This description, which addresses 

the distribution of spatial extents of thermal position- fluctuations (i.e., localization lengths) has been confirmed by 
computer simulations JtJ , and rather general, model-nonspecific arguments in favor of the broad applicability of the 
description have also been presented For any particular version of random media (e.g. the macromolecular media of 
Refs |5|,|| or the vitreous media considered in the present paper), what determines the specific content of the statistical 
description of the amorphous solid state is the form of the random constraints that the permanent covalent bonding 
imposes, and the resultant form taken by the amorphous solid order parameter. In the present context of vitreous 
media, the constraints, as we shall see below, are more intricate than they are in the macromolecular vulcanization 
context and, accordingly, the order parameter is more intricate and the statistical content more elaborate: it accounts 
not only for the heterogeneity (i.e., the distribution over the sample) of positional localization lengths but also for 
the distribution of orbital-orientation thermal fluctuations, the strengths of the position-angle thermal fluctuation 
correlations, and the statistical correlations between these physical characteristics [[)). 

In Ref. ||], in addition to introducing the model of random network forming media considered here, and formulating 
the question of the phase transition to (and structure of) the amorphous solid state via statistical-mechanical tech- 
niques, Goldbart and Zippelius made a simple variational mean-field theory for the amorphous solid state in which all 
particles shared a common localization length and all orbitals shared a common extent of their angular localization. 
The positional and angular localization parameters were then solved for, self-consistently, and it was found that, at 
a certain critical value of the density of formed bonds, a continuous transition to an amorphous solid state occurs, 
beyond which the inverse of the positional localization parameter grows continuously from zero. It was also found 
that, in response to the onset of positional localization, orientational localization of the orbitals sets in. Owing to the 
restricted form of the variational hypothesis for the order parameter adopted in Ref. |3| , specifically that it did not 
allow for the possibility that only a fraction of the particles would become localized at the transition, the critical bond 
density was over-estimated in Ref. ||. (The correct critical density was, however, known from the linear stability 
analysis of the fluid state.) Later work by Theissen et al. pQ| , in addition to allowing for networks comprising particles 
of various valencies, cured the difficulty of the critical bond density, by broadening the variational hypothesis to allow 
for a localized fraction (although it still only allowed for a single value for the positional and orientation localization 
parameters, and did not account for correlations between the thermal fluctuations of positions and orientations). 

What, then, is the nature of the amorphous solid state? If the number of permanent random covalent bonds 
introduced between particles is smaller than a certain critical value then the effect of these bonds is to bind at least 
some of the particles into random permanent molecules of a variety of types (varying in size and architecture) , each 
of which, given sufficient time, will wander ergodically through the volume of the container, i.e., the equilibrium 
state of the system is fluid. If, however, the number of bonds introduced is greater than the critical value then 
their effect is to bind a nonzero fraction of the particles into a macroscopically large disordered molecule that extends 
throughout the container, the remaining fraction of particles remaining disconnected from the macroscopic molecule 
and capable of wandering across the container, given sufficient time. By contrast, the particles that constitute the 
extensive molecule will be localized in the vicinity of random preferred spatial positions, about which their positions 
will undergo thermal fluctuations extending only over a limited spatial regime (which will vary randomly in magnitude 
from particle to particle, reflecting the random architecture of the network), and these particles will confer a rigidity 
on the entire system, so that the equilibrium state of the system will no longer be fluid and will, instead, be solid. 
Moreover, the orbitals attached to localized particles will exhibit most probable orientations, about which they will 
fluctuate thermally, the extent and nature of these fluctuations also varying randomly from orbital to orbital. In 
addition, the thermal fluctuations in the positions of the particles and the orientations of the bonds connecting them 
will be correlated, to an extent that varies randomly from particle to particle. The unconventional nature of this 
the amorphous solid state is worth emphasizing: (i) only a fraction of the particles will be localized; (ii) the mean 
positions of the localized particles will be random, as will be (iii) the spatial extent of the positional fluctuations of 
the particles, (iv) the orientational fluctuations of the orbitals, and (v) the correlations between these fluctuations 
(these parameters being characterized by a joint probability distribution); and (vi) there will be no hint of crystallinity 
beyond the shortest of lengthscales (i.e., the bond length), beyond these lengthscales the symmetries of the amorphous 
solid state being those of the liquid state. 

Our principal aims are to construct a statistical characterization of the structure and heterogeneity of the amorphous 
solid state exhibited by a model of permanently randomly bonded vitreous media in the vicinity of the solidification 
transition, and to provide a physical interpretation of this characterization. We shall do this by constructing the 
self-consistency equation for the amorphous solid order parameter, valid in the vicinity of the solidification transition, 
and obtaining an exact solution of this self-consistent equation. 

This Paper is organized as follows. In Sees. |l| and [II we shall proceed kinematically, describing the model that 
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we shall be considering, and analyzing a suitable order parameter defined in terms of the positions of the constituent 
particles and the orientations of their orbitals. Continuing kinematically, we shall explore the structure of this order 
parameter, and elucidate the physical information that it encodes. Then, in Sees. |y| and we shall address the 
model, regarding the formed bonds as quenched random variables that vary from realization to realization. By using 
equilibrium statistical mechanics, invoking the replica technique to deal with the quenched randomness, and making 
a mean- field hypothesis, we shall develop a self-consistent equation for the order parameter. By making a natural 
physical hypothesis for the form of the solution we shall 



VI, solve exactly for this order parameter in the 

near the 



in Sec. 

regime in which the thermal fluctuations of the particle positions and orbital orientations are strong (i.e. 
solidification transition). Finally, in Sec. VII we sh all ex tract from our solution a wide array of physical diagnostics 
characterizing the amorphous solid state and, in Sec. VIII we shall make some concluding remarks. We emphasize that 
throughout this work we shall be proceeding analytically, except that we shall make use of the numerically-obtained 
scaling function (of a single variable) central to the characterization of vulcanized macromolecular matter described 
in Refs. |,§. 



II. ELEMENTS OF THE MODEL 

The model of vitreous media which we shall focus on is that introduced in Ref . || , which takes as its ingredients a 
thermodynamically large number N of particles, moving in a large three-dimensional cube of volume V (on which we 
impose periodic boundary conditions) , at least some of which particles are permanently randomly bonded together 
to form a random network. At the kinematic level, the particles (labeled by j = 1, . . . , N) are characterized by their 
position vectors {cj}^^, along with the NA unit vectors {sj,a}£=ijLi describing the spatial orientations of the A 
orbitals that radiate from each of the particles j. Note that we shall be measuring lengths in units such that orbitals 
have length unity. Figure [l] illustrates the structure of the particles and the formation of a continuous random network 
out of them. 




FIG. 1. Particle structure, bond, and network formation, (a) A single particle with the near-tetrahedral equilibrium structure 
of its orbitals. In this example the number of orbitals per atom A is 4 (as would be the case, e.g., for networks of Si atoms), 
(b) Formation of a covalent bond between two particles (the participating orbitals are slightly separated, for clarity), (c) A 
collection of three particles bonded together, forming the beginnings of a random network. 

The orbitals radiating from a given particle tend to repel one another. For example, in the absence of any external 
perturbing forces, all things being equal, the orbitals of a four-orbital particle would point towards the vertices of a 
regular tetrahedron, as shown in Fig. |](a). Rather than give a detailed specification of the interactions that embody 
this orbital-orbital repulsion, we shall encode the effects of such interactions into a sequence of parameters that 
characterize the correlations between the orientations of the orbitals of a single particle. For example, we shall find 
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ourselves needing the correlator of the orientations of two distinct orbitals (<zi and (Z2) of a single particle (in the fluid 
state), say the j th : 

Y e imi ( s j^i) Y ^ni 2 ( s j,a 2 )^ i x = 7- Se u t 2 6 mi ,m 2 Ci, (2.1) 
which we have parametrized in terms of the real numbers {CV}^ (with Co = 1) that r eflect the extent to which 



the orbitals interact The angle brackets (• • which we discuss below in Sec. [IV Aj , denote thermal averaging 
with respect to a the single-particle Hamiltonian, which incorporates the intra-particle interactions. The form of this 
correlator follows from the isotropy of the distribution of the orbital orientations in the fluid state. We shall not 
find ourselves making explicit use of the correlator of the orientations of three distinct orbitals of a single particle. 
However, a simple symmetry-dictated form for it can readily be constructed, by making use of Wigner 3-j technology, 
if one wishes to compute explicitly components of the order parameter that depend on it. 

Having described the issue of a single particle and its orbitals, we now turn to the issue of the permanent random 
covalent bonds between pairs of particles, and how we are to describe them. We regard such bonding as introducing 
constraints on the relative location and relative orientation of the particles and orbitals participating in the bond. 
Specifically, we model the situation in which particles j and j' are bonded via orbitals a and a' by the constraints 
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Sj,a — Sj/ a / : (2.2b) 



7 s j,a — c j' + s j',o'j (2.2a) 



as shown in Fig. 0(b). We denote by the number M and the collection {j e , j' e ; a e , aj,}^ a specific realization of M 
bonds (i.e. a specific architecture). 

Of course, the particles in the fluid interact with one another, regardless of whether or not bonds have been 
introduced. We shall assume that pairwise interactions, depending on the relative separation and orientation of the 
orbitals, exist between all particles. The crucial consequence that we assume these interactions to have is that they 
stabilize the system with respect to the formation of macroscopically inhomogeneous or anisotropic states, such as 
regular crystalline, liquid crystalline, molecular crystalline or globular states. 



III. AMORPHOUS SOLID ORDER PARAMETER: RANDOM POSITIONAL AND ORIENTATIONAL 

LOCALIZATION 

Following the ideas of Ref. ||, which represent an elaboration of ideas discussed in Refs. fl^l , we adopt as the order 
parameter characterizing the amorphous solid state the entity 

(3.1) 

where the angle brackets (• • •) (with no subscripts) indicate a statistical- mechanical ensemble average over configu- 
rations of the particles, subject to a given collection of permanent random constraints (i.e., bonds), and the square 
brackets [• • •] indicate an average over realizations of the bonds. This order parameter, which involves products of 
replicas of a single ensemble average, depends on the collections of wave vectors {k Q }^ =1 and angular momentum 
indices {r,m Q }™ =1 . 

Let us examine this order parameter, first, in order to ascertain the nature of the physical states that it is capable 
of diagnosing, and then to understand the type of statistical information that it encodes. 



N A n 
j — 1 a—1 a—1 



A. Detection of random positional and orientational localization 



Consider the order parameter given by Eq. (3.1), and suppose that we elect to set l a = (for a = 1, . . . , n). Then 
the order parameter becomes, up to irrelevant factors of 47r: 



N n 
j=l a=l 



exp(— ik Q • Cj) 



(3.2) 
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As discussed in detail, e.g., in Ref. H, and also below in Sec. VI D, these components of the order parameter are 



capable of detecting the spontaneous random freezing of particle positions (without regard to the behavior of the 
orbital orientations). More specifically, via the wave vector dependence, these order parameter components yield 
information about the fraction of particles that are positionally localized, as well as the statistical distribution of their 
positional localization lengths. 

Suppose, instead, that we set k" = (for a = 1, . . . , n) in Eq. (jO]). Then the order parameter becomes 
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A 

-1 a—1 a 



(3.3) 



As discussed, e.g., in Refs. [| 10|, and also below in Sec. VIE2, this component of the order parameter is capable 



of detecting the spontaneous random freezing of orbital orientations (without regard to the behavior of the particle 
positions). More specifically, via its dependence on the angular indices {£ a , ro Q }™ =1 , this order parameter yields 
information about the extent and character of the orientational localization of the orbitals. It is via this compo- 
nent of the order parameter that the most direct contact is made with the Edwards- Anderson order parameter for 
Heisenberg spin glasses, which detects the random orientational freezing of magnetic moments. For example, choosing 
{t 1 , . . . , £ n } = {1, 1,0,..., 0}, and contracting appropriately on m 1 and m 2 we obtain 
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N ^ A^ 

j=l a=l 



, (3.4) 



thus recovering the familiar Edwards- Anderson form. More generally, the order parameter for random networks 
exhibits the unconventional features that the index i can be greater than unity (so that higher multipole moments of 
the distribution of orientations can be accessed) , as well as that a full characterization of the orientational freezing 
requires information from components with more than the familiar pair of thermal expectation values. 

The third category of information results from examining the components of the order parameter corresponding to 
nonzero values of both {k Q }^ =1 and {£ a }a=i- First, consider the subcase for which in every replica a at most one of 
l a and k" is nonzero. An example of such an order parameter component is 
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N ^ A^ 

j=l a=l 



^expHk 1 ■ Ci )^exp(-ik 2 ■ Cj))(^3 m 3(Sj )0 )^y/; m 4(sj 



(3.5) 



Such components measure the statistical correlations between the strengths of positional and angular localization 
across the sample. Inter alia, such components address the question: if a certain particle is strongly localized 
positionally, how likely are the attached orbitals to be strongly localized orientationally. 

Next, consider the general case in which some replicas a have nonzero values for both k Q and l a . Such components 
provide information on the extent to which positional and orientational thermal fluctuations are correlated. For 
example, as is discussed in more detail below, by setting k Q ~ in all replicas a except replicas 1 and 2, and by 
also setting i 1 = £ 2 = 1, m 1 = m 2 = 0, and £ a ~ m a — in the remaining replicas a, we would obtain access to the 
disorder average of the quantity 



(3.6) 



which is a direct measure of the extent of the above-mentioned position-angle fluctuation correlations. 

Let us pause to emphasize the three levels of randomness presented by random network forming media. There 
is thermal randomness, by which we mean the familiar thermal motion of the particles and orbitals. Then there 
is architectural randomness, resulting from the random manner in which covalent bonds are formed. Finally, there 
is micro structural randomness, i.e., the heterogeneity of the emergent solid state. This last level of randomness we 
capture statistically in a joint probability distribution that characterizes the nature of the thermal motions. 



B. Isolating the fraction of positionally localized particles 

The most basic piece of information describing the amorphous solid state concerns the value of the fraction q of 
the N particles that are localized positionally, regardless of the value of their localization lengths and the angular 
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localization of the orbitals attached to them. As shown in Refs. |3|||, this fraction q can be accessed via the order 
parameter in the following way: set l a — m a = for a = l,...,n, and then pass to the limit k — + 6 via 

a sequence for which X^a=i k° = 0- The resulting quantity is the fraction q. The reason for this is that whereas 
the value of (exp (— ik • cj)) at k = is strictly unity, the limiting value of (exp (— zk • Cj)) (as k — > 0) is unity for 
positionally localized particles, but zero for delocalized particles. For the sake of convenience, we shall refer to the 
localized fraction q as the solid fraction, and the delocalized fraction 1 — q as the liquid fraction. 



C. Distribution of positional and angular localization characteristics 



To further elucidate the physical information regarding the positional and orientational localization of the particles 
and orbi tals contained in the order parameter, we now construct a physically motivated form for the order param- 



eter (3T) in terms of certain localizational characteristics — quantities that describe the positional and orientational 



localization of particles and orbitals. 




We begin by considering the contribution from a single particle j and a single orbital a attached to it, i.e., the 
expectation value 

(e"^r/ m (s 3 , a )). (3.7) 

This function is the characteristic function of the joint thermal probability distribution describing the equilibrium 
localization of the position of the particle and the localization of the orientation of the orbital, as well as correlations 
between the thermal fluctuations between this position and orientation. First, we consider the (1 — q)N particles in 
the delocalized fraction. For such particles we have (e _lk ' Cj Y* (sj . a )) = (47r) -1 / 2 (5k, o <^,o $m.a- For the remaining qN 
localized particles we first extract from the expectation value ( |3.7|) the phase factor associated with the mean position 
fjbj = (cj) of the particle, to obtain e _ik/x j (e~* k 'fe' - '**) ^ro( s i,a))> and then express the resulting quantity in terms 
of disconnected and connected pieces: 

(e- lk - Cj Yl m ( Sj , a )) = c-' ik ^ {(e- <k -^-^)) (y/ m ( Sj> )) 

+ ^-ik-Cc-M,) _ (e -ik.( Cj -^)^ (V /m(s ., „) - (Y; m ( Sjta ))))} . (3.8) 



G 



FIG. 3. Orientational localization of orbitals. The orientation of the orbital fluctuates thermally about its most probable 
value (broken line) , the characteristic scale of these fluctuations being represented by the gray cone. 

On the right hand side of this expression, the disconnected (i.e., the first) piece contains two factors: (i) 
(exp ( — ik ■ (cj — Pj))}, which describes the positional localization of the particle (see Fig. ^); and (ii) (Yg m (sj ta )}, 
which describes the orientional localization of the orbital (see Fig. ||). If we approximate the first factor in the 
connected piece by making use of the standard cumulant expansion, by letting 3£j denote the (finite) mean square 
fluctuations / (cj — (cj)) ■ (cj — (cj)) ) in the position of the particle, and by following this strategy to second order, 
then we arrive at the approximation 

(exp ( - ik - ( Cj - fij))) « exp(-fc% 2 / 2 ). (3.9) 




FIG. 4. Orientational-positional thermal fluctuation correlations. Imagine that the particle is connected to a rather immobile 
particle by the top left orbital. As the particle moves to the shaded position its orbitals reorient accordingly. 



As for the connected (i.e., the second) term on the right-hand side of Eq. (3.8), it describes correlations between the 
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fluctuations in the particle position and the orbital orientation (see Fig. g). In the same way that we have introduced 
the diagnostic £ to characterize positional localization, we now introduce two further diagnostics: 



e- |k|2 ^ /2 rW a (k) = (( e -*-^-^) - (e- ik ^-^))(Y; m ( Sj , a ) - (Y* n (^ a )))). 



(3.10a) 
(3.10b) 



The collection of complex- valued numbers {~Si m -j,a} characterizes the orientational localization of orbital a on particle 
j; the collection of complex- valued functions {r^ m; j ja (k)} characterizes the correlations between the thermal fluctua- 
tions in the position of particle j and the orientation of orbital a attached to it. For example, consider Tio-j ja (]s). By 
expanding the exponential to first order (as we shall establish later, typical values of k are small near the transition), 
and recalling, that up to numerical factors iio(s) is s z , we see, that 



e- |k|2 ^ /2 IWk) = k • ((( Cj - - (cj - Mj)) ((b,-,.). ~ ((Hah))) + 0(k"), 



(3.11) 



which does indeed measure correlations between the positional and orientational ther mal flu ctuations, in accordance 



with the discussion of th is component of the order-parameter, given at the end of Sec. Ill A 



By rewriting Eq. (3.8) in terms of these diagnostics, and making use of the approximation ( |3.9| ), we arrive at the 
form 



Yf m (.H°)) ~ e ~ ik ' Mj e-i k i 2 « 2 / 2 {s, ma , a + r /m!j|0 (k)}. 



(3.12) 



By inserting this form, appropriate for particles j that comprise the localized fraction, into Eq. (3J), and incorporating 
the contribution from the delocalized fraction, we arrive at the form 



^ N ^ A n 

n 

j — 1 a—1 a — 1 



-ik" 



.(Sj,a) 



(i-?)n^k-,o 



1 , A n 

EiEn^^- |k " |2c?/2 {^ 



N ^ A 

j loc. a—1 a—1 



•>(k a )} 



(3.13) 



= (l-q)f[6 k a, <Q + q J d 3 fi J dTd{£}V{T}P(»,T,{X},{T}) 



Tl 

;e -^-E: =1 ^ e -^E: =1 i ka i 2 



a(k Q )}, 



(3.14) 



where we have made the definition r = l/£ 2 . The integration measures and £>{r} respectively denote the 

multiple measure W.t m dHi m and the multiple functional measure Yiem'^ > ^ im - We nave a l so introduced the joint 
probability distribution P, central to our characterization of the localized particles in amorphous solid state, defined 
via 



p(M,r,{£},{r}) 



oo t 
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TV ^— ' A 

j a=l 



fcOra=-< 



(3.15) 



in which the final factor D(-) is a functional delta- function. 

The next step in our construction of a physically motivated form for the order parameter involves the identification of 
specific symmetries that we anticipate the amorphous solid state to possess, viz., macroscopic translational invariance 
(MTI) and macroscopic rotational invariance (MRI). MTI reflects the notion that although in the amorphous solid 
state translational invariance is spontaneously broken at the microscopic level (in any particular realization of the 
disorder) , this invariance is restored at the macroscopic level, in the sense that no quantity computed by averaging over 
any macroscopic subvolume of the system exhibits any dependence on the particular subvolume chosen. Similarly, MRI 
reflects the notion that although rotational invariance is spontaneously broken at the microscopic level, it is restored 
at the macroscopic level in the sense that no quantity computed by averaging over any macroscopic subvolume of the 
system exhibits any orientational preference. 
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As for MTI, it amounts to the hypotheses: (i) that the disorder-averaged distribution P exhibits no correlation 
between the mean location of a particle and its other statistical characteristics; and (ii) that the distribution is 
translationally invariant (i.e., is independent of fi). 

Although MM also imposes conditions on the joint probability distribution P we do not need to impose these 
conditions explicitly. The reason for that is that MM for the component i7(k; 0, 0) is assured by the fact that £oo;j,a 
and roo ; j,a are constants (in fact, one is zero), and the assumption that the localization clouds of the particles are 
spherical and, accordingly, described by the single r.m.s. value of the fluctuation in the particle's position 3£|. MM 

for the anisotropic components of f2(k; £, m) is a consequence of the fact that they, as we shall see below, are perturbed 
away from their zero values by MM-satisfying couplings to the 0(k; 0, 0) component. Thus MM is assured by the 
theory itself, and does not need to be explicitly incorporated into the proposed form of the order parameter. 
By making use of the MTI hypotheses we arrive at the form: 

P(M ; r,{S},{r»= P(r ' { ^' {r}) , (3.16) 



which, when inserted into Eq. (3.14), leads to the expression 

n roc p 

(4^)"/ 2 fi(k; I m) = (1 - q) TT <5 k «.o &-,o <W,o + g V k- o / dr \ d ^ V ^ P ( r ' { S >> { r » 

a— 1 J 

n 

xe -*SI = oi kQ i 2 JT {s^ mQy - a + r, QmQ;J - Q (k Q )}, (3.17) 



a=0 



where we have introduced the notation f2(k; £, m)| k o =0 ^o =m o = o to denote the order parameter ( |3.1| ), and where hats 
indicate n+ 1-fold replicated versions of quantities. By extending the result of this approach to include the dependence 
on the zero-replica variables (k°,£°,77i ) in a permutation-symmetry-dictated way we arrive at the form 

(4^)( ,l + 1 )/ 2 n(k;U) = (1 - 9 ) p 4-,o ^«,o^«,o + ?i i|0 1 dr J V{T}p(t , {£}, {r}) cr^ 2 !^ 

n 

U {^«m«;i,ii + ^e a m a ;j,aO s °')}- (3.18) 

As we shall see below, for our solution we shall need an assumption for the form of the order-parameter component 
fl(k; 0, 0). To motivate this assumption, we set I — rh = in Eq. (3.18). As is easy to see from the definitions (3.10a) 



Q=0 



X 

a=0 



and (3.10b), Soo y ,a = l/v47T and Too,j,a = 0, leaving us with the form 

71 />oo 

(4^)( n+1 )/ 2 fi(k;6,6) = (l- ( ?)n^,o + ^ ko / dT P (r)e- k2 ^ 7 (3.19) 



where p(r) is a reduced form of the full joint probability distribution P(r, {£}, {T}), and describes only the positional 
localization of the particles: 

p(r) = J V{T} P(r, {E}, {r». (3.20) 

We note that, up to trivial factors of An, the expression ( 3.19 ) is identical to the Ansatz used in Refs. [|]||] in the 
context of vulcanized macromolecular media. 



IV. DISORDER-AVERAGING; REPLICA STATISTICAL MECHANICS 

Having described the relevant "kinematics" i.e., the degrees of freedom and the constraints that characterize the 
model of randomly covalently bonded particles, we now formulate the statistical mechanics of the system, paying 
particular attention to the quenched (i.e., nonequilibrating) nature of the random constraints. At this stage we shall 
be following the method sketched in Ref. [|3| which itself builds upon the general approach to macromolecular networks 
introduced in Ref. 
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A. Partition function 



The partition function of the system, subject to the constraints {j e , j' e ] o, e , a' e } e l , which we collectively denote by 
C, relative to the partition function of the unconstrained system, is given by 

M 

Z(C) = ( [] 5 (3) (c, c + 2*»J.,». - - ~ Sj , K ) A^(s je<ae + s jLl< )) N x . (4.1) 

e— 1 

The angle brackets denote equilibrium averaging with respect to a Hamiltonian that incorporates interactions between 
distinct particles, as well as between the orbitals of a single particle. The subscript indicates that this average is taken 
over one copy of a system of N particles, and anticipates the introduction of replicas that we shall make shortly. The 
two types of delta-function, 5^ and A^- 2 \ serve to eliminate configurations that fail to satisfy the constraints, and 
are appropriately defined in the following way: 

*(3) (C1 - c 2) = ^) (C2 -co=s {^^) {^^y> ^ 

oc I 

A( 2 )( Sl ,s 2 ) = A( 2 )(s 2 , Sl ) = £ Yl ^m(si) Y/ m (s 2 ), (4.2b) 

£=0 m=-l 



where, (corresponding to the periodic boundary conditions imposed on the system), the sum over k is taken over the 
Cartesian components k v — 2T:n u /V 1 ^ , with n v being integers (both positive and negative), and the Ye m are the 
usual spherical harmonic functions, the arguments of which are unit vectors. As one can readily check by making use 
of the orthonormality properties of the spherical harmonic functions, i.e., 

3 /expiki • c\ * /expik 2 • c\ 
d c Try Ttt = S *ite> ( 4 - 3a ) 



v ) \ Vv 

d2 sY; imi (s)Y e2rri2 (s) = 8 tui2 5 muma , (4.3b) 



A' 



where the subscript V on the former integral indicates integration over the volume of the sample and the subscript S 
on the latter integral indicates that the integration is taken over the two-dimensional surface of a three-dimensional 
sphere of unit radius, these delta-functions possess the basic properties: 

d 3 c 2 <5 (3) ( Cl - c 2 ) 5^ (c 2 - c 3 ) = <5 (3) ( Cl - c 3 ), (4.4a) 

/ d 2 S2 A( 2 >( Sl ,s 2 )A( 2 )(s 2 ,s 3 ) = A (2) ( Sl ,s 3 ). (4.4b) 

Strictly speaking, the partition function Z is correct only up to the appropriate Gibbs factorial factor, and would 
not, as it stands, give rise to an extensive free energy hence the tilde. As we shall be focusing on the order-parameter 
self-consistency equation, in which (as is well known) the Gibbs factor plays no role it can be safely omitted here. For 
a detailed discussion of this issue, see Sec. 2.4 of Ref. M. 



B. Deam-Edwards distribution 



At this stage, we introduce a statistical distribution characterizing the realization of the random bonds, following 
the elegant strategy of Deam and Edwards jl3) . We take for the probability density that the collection of bonds C is 
formed the quantity 

Vm\C) cx — Z{C), (4.5) 

which is analogous to the probability density chosen by Deam and Edwards for the case of vulcanized macromolecular 
networks [ p^[ . Instead of working with a fixed number of bonds, we allow their number to vary in a quasi-Poissonian 
way, and control the mean number of bonds by the control parameter /i 2 . For a given value of M, the Deam-Edwards 
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distribution is proportional to the probability density for finding the set of pairs of orbitals {jm> 3m'i Q>'m 
be overlapping. The factor /j, 2M represents the probability that a bond is formed between each of these M overlapping 
pairs. Thus, the Deam-Edwards distribution provides a statistical characterization of a process of forming permanent 
bonds in which constraints are instantaneously introduced into the liquid state at equilibrium. As such, it is an 
idealization of the random-network-forming process, which generally takes place on a time-scale during which at least 
some relaxation of the structure can occur. To handle the complication of relaxation would require the introduction 
of kinetics into the description, rather than purely equilibrium notions. Said another way, one can view the Deam- 
Edwards distribution as a strategy for freezing in liquid-state correlations, as process that is regarded as happening 
spontaneously in glass-forming systems, but here is introduced externally. The distribution encodes the physically 
attractive feature that the networks it gives appreciable weight to exhibit the macroscopic symmetries of the liquid 
state, inasmuch as the bond collections to which it gives appreciable weight correspond to likely configurations of the 
liquid state. With this distribution of bonds some of the correlations of the liquid state are quenched in, to a degree 
controlled by the mean number of bonds formed. 



C. Replica representation of the disorder-averaged free energy 

We now set about constructing the disorder-averaged free energy per particle (relative to that of the system prior 
to random covalent bonding) /, which is defined via 

Npf= [lnZ(C)], (4.6) 

where j3{= l/k^T) measures the inverse temperature. As mentioned in the previous subsection, the Gibbs factor has 
been omitted, but this will have no consequences for the order-parameter self-consistent equation. By making use of 
the replica technique (see, e.g., Ref. @), / can be obtained via: 

/=lim/„, (4.7a) 

n— >0 

-npNf n = [Z n ] - 1 = {Z n+1 - Z x ) jZ x , (4.7b) 



Zn+1 S ( 1 E Ell ^(C? + ¥la ~ $ ~ ¥1>,a>) * i2) Ka, -Sji,,) ) ) (4.7C) 

\ V j,j'=l o,o'=l a=0 / / N,n+1 

Here, Z n+X i s th e replicated partition function, arising from the averaging of Z n over the Deam-Edwards-type 
distribution ( |4.5| ), and the denominator Z x arises from the normalization of the Deam-Edwards distribution (see 
App. ^ for details). Notice the striking occurence of a theory involving n + 1, rather than the usual n, replicas, a 
feature, arising from the presence in the partition function in the Deam-Edwards distribution. (The extra replica 
"computes" the distribution of quenched random bonds.) The angle brackets (• • -) n +i indicate an n+ 1-fold replicated 
normalized average over the positions of the particles and the orientations of the orbitals, weighted suitably by a 
Hamiltonian that does not couple the replicas. 



JV 



As one can see from the exponent in Eq. (4.7c), the replicated theory possesses the symmetries of independent 
translations and rotations of the replicas, i.e., 

c" -^TZ a -c a + a Q , (4.8a) 
s2-ft°-s» (4.8b) 

where {a"} are n + 1 independent arbitrary translation 3- vectors, and {lZ a } are n + 1 independent arbitrary 3- 
by-3 rotation matrices. As we shall see, the transition to the amorphous solid state is marked by the spontaneous 
breaking of the symmetries of the relative translations and rotations of the replicas; the common translations and 
rotations remain as residual symmetries. These residual symmetrie s corr espond to the macroscopic translational and 
rotational symmetry of the amorphous solid state discussed in Sec. [II C| . The theory also possesses the symmetry of 



the permutation of the n + 1 replicas, which remains intact in the amorphous solid state. 
For the sake of convenience, we introduce the replicated delta-functions, defined by 



fax-ca^lT d&(^ -eg), (4.9a) 
A(§i, s 2 ) = JJ™ A( 2 )K,s«), (4.9b) 
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where c denotes {c°, c 1 , . . . , c"} and s denotes {s°, s 1 , . . . , s"}, and also the replicated spherical harmonics Y, defined 
by 

y^(fi) = U n ( 4 - 10 ) 

where £ and m respectively denote {£°, i 1 , . . . , £ n } and {m°, m 1 , . . . , m"}. 

V. MEAN-FIELD APPROXIMATION 
A. Self-consistency condition for the order parameter 



We now develop a mean-field approximation for the replica partition function, Eq. (4.7c). To do this, we rewrite 
the partition function as follows: 



N A 



Z n +i = ^exp \2ttNVh 2 J^dc jf dl ^ ^ £j*((cj + Js i>0 ) - 6) A(S ilOJ s) 

1 N A v\ 
X ^EE + 3 8 i'.»0) A(s,-%',a') J ) (5-la) 

= ^cxp (2ttNVh 2 J v dc dsj^j Yl E - ( £ J + 2 § ^)) A ( § ' 

1 N A s\ 

x ]^E £*( £ -( £ i' + ^«0)A(-M^a')) ) , (5.1b) 

i'=la'=l ' iV,n+l 

where J v dc denotes n«=o Iv d 3 c a , and J s ds denote s TT!!-n fg d 2 s a , and where, to obtain the last form, we have used 
the symmetry property of the delta- functions, Eqs. ( 4.4a , 4i4q ), 



Next, we introduce the (real-space version of the) amorphous solid order parameter, 

I n A y 
m §) ^ ( NA E E 5 ( £ - & + K'* ! ) A ( § , Si,a) ) ■ (5-2) 

\ j=l 0=1 / 

Then, upon setting 

^ at A 

am E E - ( £ j + 2 § j.«)) A ( § ' = n ( £ ; § ) + m ^ § )' (5- 3 ) 

i=l a=l 

i.e., the order parameter 0(c;s) plus the fluctuation #0(6; s), expanding the exponent in powers of <50(6;s), and 
omitting terms quadratic in $0(6; s), we obtain 

Z n +i ~ exp ^— 2nNVy? J dc J ds 0(6; s) 0(6; — s) 

-h/Vln^cxp ^AirVn 2 J dc J ds 0(6; -s)^ ^ 5(6 - (6 X + £ii )0 )) A(b, s m )^ J , 



(5.4) 



where the resulting expectation value involves only the position and orbital-orientations of a single particle. 

Th e re ader will have observed that the mean-field approximation strategy has yielded a one-particle problem, 



Eq. (5.4), as desired. However, there is a subtlety associated with the manner in which the various interactions 



present in Eq. (5.1t) are treated, which we now address. The angle brackets in Eq. (5.1b) denote averaging over n + 1 
(coupled) replicas of the N (coupled) particle system. The intra-replica coupling originates in the interactions between 
particles present in the liquid state; on the other hand, the inter-replica coupling originates in the random constraints. 
As discussed in detail in Sec. 5.1 of Ref. it is useful to transfer the so-called one-replica sector contribution to 
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the inter-replica coupling to the intra-replica coupling (which is thereby renormalized) . (The intra- and inter-replica 
couplings both contain trivial contributions in the zero-replica sector, as does the order parameter; we ignore these 
contributions.) The subtlety is that the structure of the theory in the one-replica sector is quite different from that 
in the higher-replica sectors: whereas the constraints tend to destabilize all sectors, this tendency is counteracted in 
only the one-replica sector by the original intra-replica interactions. Consequently, at the amorphous solidification 
transition the one-replica sector component of the order parameter remains zero, whilst the higher-replica sector 
components become nonzero. Indeed, the competition between these two processes can be r egard ed as a form of 
frustration, which resolves itself by the formation of a state possessing MTI and MRI (see Sec. [II C ). 

On a technic al level, this discussion amounts to the following dictum: in all subsequent equations, e.g., 
Eqs. (|5.4 5.5,5.7), the component of the order parameter lying in the one-replica sector is to be set to zero. Ac- 



cordingly, the self-consistent equations that follow pertain to all sectors except the zero and one replica sectors. This 
notion is straightforward when the order parameter is expressed in the (plane and spherical) harmonic representation 
[as it is, e.g., in Eq. fl5.7| )]. In this representation, setting the one- replica sector contribution to zero refers to setting 
to zero the contribution in which nonzero entries in {k; £, m} appear in precisely one replica. (By the zero-replica 
sector we mean the sector with k = £ — m — 0.) 

We now return to the task of obtaining a self-consistent equation for the order parameter. By making the partition 
function (5.4) stationary with respect to f2(c;s) we arrive at self-consistent equation (SCE) for the order-parameter: 



' l A 

— ^2 K^o - (ci + £si, tt )) A(s Q , Si, a ) 



x exp ( 4ttV> 2 / dc [ ds fl(c;-s)i Vi(c- (ci + pi, a )) A(s,s M ) ] ) 

V JV JS 2^=1 / / l,n+l 

exp ^4ttV> 2 J^dc J ds fi(c; -b)~ £j §{c - (c x + ^s li0 )) A(s, s x , ) J J 



(5.5) 



The presence of A-function factors provides the option of replacing the dynamical variable si )0 in the argument of 
the ^-functions by the parametric variable s, which replacement we sometimes make. 

The transformation of the order parameter to a representation in terms of the plane wave and spherical harmonic 
coordinates, and the inverse transformation, are effected as follows: 



il(k;£, m) = / ds (s) exp(Aik • s) / dc exp(— zk ■ c) f2(c; s), 
■Is "" " Jv 

§ ) = J2i ex P( zk • £ ) exp(-|ik • s) Y m(*) n(k; £, m) 



(5.6a) 
(5.6b) 



This choice of transformation has the effect of keeping the physical interpretation of the order parameter f2(k;£, m) 
free of any complicating factors of exp (^ik • s), leaving the full factor exp (ik ■ s) in the fluctuating variable to which 
the order parameter couples. Via this transformation, one arrives at a transformed self-consistent equation for the 
order parameter: 



tt(k ;£ ,m ) 



47171 



Eti e^ k °- c y? Ao (S B ) exp ( ^ Y: U .rn ^ ^ ™) e *- C 3 E:=i e lk ' s « Y^(-s a ) 



l,n+l 



exp 



1,71+1 



(5.7) 



We shall often refer to the triple (ko;£ ,mo) as tne "external variables" 



B. Instability of the fluid state 



We now demonstrate that upon increasing the density of formed bonds the fluid state is rendered linearly unstable. 
To do this we follow Ref. ]l0| , and expand the replica free energy, n/3Nf n in Eq. ( 4.7b| ), to second order in the order 
parameter Q, thus obtaining 
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k £1 ,mi ^2)7^2 

^47r/i 2 y ivy 

It i rh CLl ,a2 a — 



) 2 |l e lofk;*.*.)! 2 (-if e n ( *-■ +(i 4 -;-- )c " ). (5. 8 ) 



where, as discussed in Sec. 0, the coefficients Ci represent the effects of the interactions between the orbitals of a 
single particle (and are all less than 1 for I > 1), and the kernel M is given by 



(imi ,1-zrhi * 



,£)=/ rfS ( § ) %h 2 (- § ) CX P(^ ' § )- (5-9) 
Js 

We remind the reader that this linea r stability analysis applies only to the higher replica sectors of the order parameter, 



as per the discussi on in Sec. [V A| . The corresponding analysis applied to the one-replica sector reveals the fact, 
anticipated in Sec. VA that the one-relica sector remains stable at the transition, owing to the stabilizing effect of 
the inter-particle interactions. 

We expect the liquid state to become unstable first for long wavelengths, corresponding to k — > 0. In this limit, 
M| iA | „ (k) — > | 2 £mi,m 2 ( — l) e+m . By examining the coefficients of the quadratic terms in Eq. fl5.8j), and specif- 
ically their signs, we see that for fi 2 < /.i 2 = 1 the coefficients are positive for all components of the order parameter 
[i.e., for all values of (k, £, rh)\ and, therefore, that the free energy has a local minimum at — 0. Thus, for /i 2 < fi 2 
the fluid state is (at least locally) thermodynamically stable. On the other hand, for fi 2 > /z 2 , certain coefficients be- 
come negative, starting with the longest length-scale (and isotropic, corresponding to I = 0) modes. This sign-change 
indicates the loss of the (linear) stability of the fluid, and the concomitant acquisition of a nonzero value of the order 
parameter. As usual, the linear instability of one state does not sharply specify the nature of the stable state that 
replaces it, although the directions of instability do provide hints. In the present setting, the residual stability of the 
one-replica sector suggests that the primary characteristic of the new state is macroscopic translational and rotational 
invariance, which is the mechanism by which the induction of energetically-costly order in the one-replica sector is 
avoided. Thus, it is reasonable to anticipate that the state that replaces the fluid state upon the formation of a suffi- 
ciently large density of bonds is the amorphous solid state. Furthermore, as the coefficients Ci are smaller than unity 
for I > 1, and become progressively smaller with increasing £, we may conclude that all anisotropic sectors remain 
stable, at least for bond densities in the vicinity of the amorphous solidification transition. Thus, it is reasonable to 
anticipate that anisotropic ordering (i.e., orientational localization) will arise only as a response to positional ordering, 
via nonlinear coupling between isotropic and anisotropic order-parameter components. Thus, as will be borne out 
below, we should anticipate that the form and extent of the anisotropic ordering will be computable algorithmically, 
as a perturbative correction to the nonperturbative result for the ordering in the isotropic sector p3 . 

VI. SOLUTION OF THE ORDER-PARAMETER SELF-CONSISTENT EQUATION 

A. General strategy 



Having obtained the self-consistency equation (SCE) for the order parameter, Eq. (5/7), we now turn to the issue of 
solving it. We shall begin by extracting from the full SCE a transcendental self-consistency equation for the fraction 
of localized particles q by considering the SCE at external variables £q = rho — and taking the limit ko — > 0. Solving 
this equation in the vicinity of the amorphous solidification transition (i.e., for small excess crosslink densities) we 
will find that q tends to near the transition, allowing us to expand the SCE for the order-parameter in powers of q 
and truncate the expansion, retaining terms of order q 2 . 
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FIG. 5. Component- and band-structure of the order parameter (see text for explanation). 

We will then solve the SCE for individual order parameter components, starting with the only component with an 
unstable band, the isotropic component f2(k;0,0) (the unstable band being those long wavelength modes for which 
2e), which we shall obtain by utilizing the form of the solution, due to refs. (|||, and discussed in Sec. Ill C 



The remaining (anisotropic) components are linearly stable and thus we can solve for their leading-order values 
perturbatively, by considering their couplings to f2(k;0,0). We shall obtain these anisotropic components by first 
solving for the leading-order contribution to the two lowest angular-momentum components, and then obtaining the 
solution in the general case by induction. The structure of the order parameter is illustrated in Fig. ^[ showing the 
different components and bands. 

Throughout the entire calculation we shall only be concerned with finding the leading-order contributions to the 
components of f2. This will generally imply (i) truncating the expansion in powers of q (typically at second order), 

(ii) ignoring the coupling of the components of the order parameter to higher angular-momentum components, and 

(iii) truncating expansions in powers of k (typically at linear order), as we expect typical values of k 2 to be of order 
e. Many technical details of the calculations have been relegated to the appendices. 



B. Fraction of positionally localized particles 



The first step in our soluti on of the order-parameter SCE is to determine the fraction of localized particles q. 
Following the discussion in Sec. Ill B , we first separate the delocalized and localized fractions in the full order parameter 
fl by writing 



ft(k;f»E(l- 9 )nf^ 



a=0 



,0 »m»,0 



4tt 



qW(k;e,m), 



(6.1) 



where qW(k;£, to) is the part of the order parameter describing the localized particles and is analytic at the origin, 
with W(0; 0, 0) = I/vItt The delocalized contribution cancels from the numerator and denominator of the SCE, so we 
may rewrite Eq. (5.7), replacing fi(k;^, to) with qW(h;£, to) in both the numerator and the denominator, obtaining 



i E ^ =i ^ko-c (s Q ) exp £u* q W(k; I, to) -\ ^Li f ^-^)) 1<n+1 

\ ' ' / / l.n+1 



(6.2) 



Following the ideas of Refs. HH, to obtain a SCE for q we consider the SCE (6.2) for the order-parameter component 
ft(k ; 6, 6) in the limit k 6. We start with the SCE for fi(k ; 6, 6): 



O(k ;6,6) =c" 



fj, q 



1 71+1 / A 2 1 



■in 



(6.3) 



l,n+l 
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Note that the denominator of the right-hand side of Eq. ( |6.2| ) has been replaced by exp(/i 2 g) (see App. |c] for details). 
We now consider the limit k — ► via a sequence for which k = 0. The left-hand side, as can be easily seen from 
Eq. ( |6.1[ ), becomes qj\pkn. To evaluate the right-hand side we follow the procedure used in App. ^| to obtain the 
denominator of the right-hand side of the SCE. 

We expand the exponential in a power series in q and consider the r-th order term, and pass to the replica limit 
n -> 0: 




... W(ki;4,mi) ...W{k r ;£ r ,m r )e lk ° c e lklC ...e ik *" c 

U. r £ r m r 

J2 e ikl ' s » 1 ...e ik ^Y &mi (-s ai )...Y frmr (-s ar )\ . (6.4) 

,...,Or I l x 

By MTI of W(k; £, rh) we know that k^ = 0, which means that kj = 0. Also, using MRI of W (k; rh) we see that if 

k = 6 then ( |C3| ) can only be nonzero if £ = 6 also. Knowing that W(6; 6, 6) = for the r-th order contribution 

(with r > 1) we get — ] {n 2 ) r q r . For r = 0, however, we get 0, becausethe limit of (exp(— iko • c)) as ko — > is zero. 
Resumming the power scries we finally obtain 

( e ik °' e exp 1 W & l > ™) e *' 6 T E ^''"^aC 8 -)) ) = E ^ 0* 2 )V = ex P^ 2 «) ^ L ( 6 - 5 ) 



km a=1 ' i - ' i r>1 



n 

1,71+1 



We thus find the following self-consistency equation for the fraction of localized atoms: 

l-q = cxp(-/i 2 (j). (6.6) 
This self-consistency condition is precisely that obtained in the case of vulcanized macromolecules [ppl, and earlier, 



in the context of random graph theory, by Erdos and Renyi (see Ref. |6J). 

We will find it convenient to exchanging the control parameter /i 2 for e, defined via /i 2 = 1 + e/3, so that e vanishes 
as the transition is approached. The self-consistent equation for q is transcendental, but it is easy to analyze it 
graphically. Then, for e < we find that there is only one solution q = 0, corresponding to the liquid state of the 
system (no localized particles), and for e > we find that for e small q is small also, indicating that the fraction q of 
particles comprising the amorphous solid state tends to zero in the vicin ity of the amorphous solidification transition. 



We can thus expand the exponential on the right-hand side of Eq. (6.6) obtaining the fraction q to first order in e: 



g«-e. (6.7) 



C. Perturbation expansion for the self-consistency equation 



Having found the fraction q of localized particles to be s mall in the vicinity of the transition, we may expand the 
self-consistency equation for the order parameter, Eq. (6.2), in powers of q, to second order, obtaining: 



qW(k ;£ ,m ) + {l-q) J] 



<5k a ,0 <^<*,0 ^m a ,0 



4tt 



-fj, q 



+ 



.4 



E e ~' ko - 6 ^Uo^) 



ao = l 
A 



\A 

e -M 2 ? / 1 



ao = l 
A 



1,71 + 1 

2 



ki ,£i ,771 1 



ai=l 
A 



l,n+l 



a ° = 1 ki,£i,rfii a i = 1 



(6.8) 



1,71 + 1 



16 



Cancelling the liquid contribution to on the left-hand side with the 0-th order contribution on the right-hand side, 
rearranging terms and replacing exp(— /J, 2 q) with 1 — q we arrive at the form of the SCE that we shall focus on: 

A 2 1 ^ 

?w(Mo,rf»o)«(i E wfahM^ E (^' e ^( s ».) sitl ' 8(iitl4 ' i 4(^)) w 

kiJi.mi a ,ai=l 



1 — q f Att/J, 



2X 2 



2 V ^ 



n 



I 2 



E E ^(ki;4,mi)iy(k 2 ;^ 



ki,<i,mi k 2 ,^2,™2 



x 73 E ( e "*°" e* 1 ' 6 e^- § - Y lrh {-s ai )^ ^(-Sa,)) ■ (6.9) 



ao,ai,a2 — 1 



D. Self-consistency equation: Isotropic sector 

Having obtained the localized fraction q and verified the consistency of the expansion of the SCE, we now turn to 
the issue of solving th e SC E for individual components of the order-parameter. 



As discussed in Sec. VB , near to the amorphous solidification transition the only linearly unstable band of W is the 
long- wavelength band of the isotropic component W(k) (i.e., W(k; £, m)|| =7fl _g having sufficiently small (replicated) 
wave vectors, specifically, k 2 < 2e). Thus, the basic process occurring at the transition is the acquisition of a nonzero 
value by the unstable components of the order parameter, which in turn perturb the stable components away from their 
zero values. The stable components include both the band of the isotropic component for k 2 > 2e and the anisotropic 
components for all values of k. As our aim is to calculate the leading contributions to each of the components of 
W at small positive e, we may, as a first step, proceed by computing the self-consistent value of W(k), neglecting 
the feedback on its value coming from the nonzero values that it induces in the (stable) anisotropic components 
[i.e., W(k; £, rh) ||^g] . Although the k 2 > 2e band of the isotropic component W(k) is linearly stable, it is necessary 
to treat it self-consistently together with the unstable band, owing to the fact that they constitute a continuum (see 
Fig. H), and therefore the stable band includes elements with an arbitrarily large "susceptibility" to perturbations 
caused by their couplings to the elements of the linearly unstable band [i.e., W(k) for k 2 < 2e\. 

We therefore consider the SCE for W(k\ £,m), Eq. ( |6.9D , for isotropic external arguments (i.e., £ = rh = 6), and 
ignore the effects of all anisotropic components on the right hand side. We thus arrive at the closed, nonlinear SCE 
for W(ko): 



W(k ) = \l - | - L j W/(k ) + |^E kl ^(ki)W(ko - kr). (6.10) 

Precisely this equation emerges in the context of randomly crosslinked macromolecular networks, from both semi- 
microscopic and Landau-type approaches |^,^[. In that context, the order parameter has only isotropic components, 
in contrast with the present context. To solve Eq. ( |6.10D we invoke the hypothesis for W(ko) discussed in Sec. Ill C , 



viz., a parametrization in terms of the normalized probability distribution of localization lengths p(r), along with the 
factor, enforcing MTI of the solution: 

W(ko) = -^= rfT P (r)e- k! / 2T . (6.11) 
V47T J a 

As is shown in Ref. ||, this leads to the following nonlinear integro-differential SCE for p(r): 

i.e., precisely the equation found in Refs. By making the rescalings 

r -> 9 = 2r/e, (6.13a) 
p(t) -> tt(0) = |p(r), (6.13b) 
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as in Refs. 
equation 



[||]6), we determine that the universal scaling function n(6) satisfies the nonlinear integro-differential 



Yd9 



(1-0)tt(0) 



dd'iri 



(6.14) 



together with the normalization condition 
numerically, and the result is given in Rcfs. 
results from numerical simulations in Rcf. 



jT 3 d9n(9) = 1. The resulting scaled distribution ir(9) can be obtained 
|||. (The prediction for this universal scaling function is compared with 



.) Thus we have obtained the isotropic component iy(k) of the order 
parameter W(k;£, rh) in the vicinity of the transition. As we have discussed in Sec. Ill A , we are thus in possession 
of statistical information concerning the spatial localization of particles, regardless of the angular localization of the 
orbitals. 



E. Self-consistency equation: Anisotropic sectors 

We now turn to the task of calculating the leading-order contributions to the anisotropic components of IU(k; i, rh) 
[i.e., W^(k; £, to)||jq]- We remind the reader that the anisotropic components are all linearly stable near the transition, 
and not merely infinitcsimally so (i.e., linear stability analysis indicates that none of these anisotropic components even 
become marginally stable at the transition). Note the contrast with the stable band of the isotropic component |W(k) 
for k 2 > 2e], which, though stable, do include marginally stable components (i.e., components of arbitrarily small 
"mass" ) . Unlike the stable band of the isotropic component the (also stable) anisotropic components are separated 
by a "gap" from the unstable components, owing to the discreteness of the external variable £, the components of 
which take on integer values only, in contrast with the components of k, which are continuous (in the thermodynamic 
limit). It is this "gap" that allows us to obtain the stable anisotropic components by means of perturbation theory, 
which we could not use to solve for the stable band of the isotropic component (see Fig. ||). 



1. Anisotropic sector: Angular momentum 1 in one replica channel 

Rather than begin with generalities, we first consider the lowest angular momentum sector {^(k; £, m)}^ 2=2 (i- e -' 
the collection of order-parameter components W(k;£, rh) such that £ 2 = '}2 a £ a {£ a + 1) = 2. In this case £ is some 
permutation of the form (0, 1, 0, 0, . . . , 0). We therefore consider the SCE for W"* (k), i.e., Eq. (|6.9|) for the anisotropic 
external argument 



replica 1 ... ot\ a\ + 1 ... n 

£ /0, 0, 1, 0, 

rh I 0, 0, . . . , mi, 0, . . . , 



and arbitrary k, which, as we establish in App. |^, reads 

WT mi (k) = C (0) (-WZ, (k) + ^(ki) C t * - ^ W ^ ~ k i) W i™ ( k i) 



(6.15) 



(6.16) 



The symbol k^* denotes the complex conjugate of the m-th spherical tensor component of the vector k a spherical 
component of the vector k a [see App. Eq. (A3), for the definition]. The parameter encodes physical information 
arising from orbital-orbital correlations of a single particle and is defined in App. [d]. Specifically, depends on the 
£ = 1 value of the free-particle two-orbital orientation correlator (Y e * m (s a ) Yii m i (So/)) x,x- The permutation symmetry 
among the A orbitals and rotational invariance of the joint probability distribution of their orientations restrict this 
correlator to have the form 

(y/ m ( So )^ W (s o 0) ii = ^<VSm,m'(<W + (l-<W)Q), (6.17) 

characterized by the parameters {Ci\fL x (with Co = 1, b y normalization), introduced in Sec. |l} As for the issue 
of what terms have been omitted in arriving at Eq. J6.16D , we are concerned only with the leading-order values of 
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the components of W, and therefore, here and elsewhere, shall omit terms that do not alter leading-order values. In 
particular, in arriving at Eq. (3.16) from Eq. (|S.S|) we have neglected all components of W having £ 2 > 2. Such terms 



are sufficiently small that feedback from them would not alter the leading-order value of W"^ (k) . We shall verify the 
internal consistency of this assumption below. In addition, we have only kept terms carrying sufficiently few powers 
of components of k. A s th e charact erist ic value of k 2 is of order e, higher powers render terms subdominant. The 
steps leading fromEq. ( |6.9[ ) to Eq. ( |6.16| ) are explained in detail in App. |j. 

To solve Eq. ( 6.16 ) we rewrite it in the form of a (Type II inhomogeneous) integral equation by moving the first 
term on the right hand side to the left, regarding the second term as a known inhomogeneity, and the third term as 
a perturbation. Wc then solve this equation iteratively, thus obtaining 



w^ i (k) = -J=c« w(k)k«\*, 



(6.18) 



where the parameter (JW depends on A and C\ (see App. [5]). In fact, for this particular component of W our 
procedure merely amounts to truncating the Born series after the zeroth order (i.e., effectively ignoring the perturbation 
altogether), and solving the resulting algebraic equation. However, for certain other components it turns out to be 
necessary to retain the first-ord er ter m, for reasons that we shall explain below. 

As we have discussed in Sec. [II A , the result that we have just obtained about the value of the order-parameter 
component (k) yields statistical information concerning the variations, across the system, of the strength of the 

correlations between the thermal fluctuations of the positions of the localized particles an d the t hermal fluctuations 
of the orientations of their orbitals. We will discuss this information in more detail in Sec. VII C 



2. Anisotropic sector: Angular momentum 1 in two replica channels 

Having obtained the leading-order contributions to the order-parameter components W(k) and W"^ (k) [i.e., the 
isotropic (largest) and anisotropic (next largest) components], we now address the component of W corresponding to 



replica 1 ... ol\ ol\ + 1 ... n 

£ /0, 0, 1, 0, 

m I 0, 0, . . . , mi, 0, . . . , 



(6.19) 



which we denote by W"^™^ (k). The motivation for examining this component is that, in contrast to the components 
W(k) and (k), the limit k — * of this component provides information purely about the orientational localization 



of the orbitals, independent of the positional localization properties, as we shall discuss in more detail in Sec. VII C 

As shown in App. [^, the retention of all terms that give rise to leading-order contributions to (k) leads 

to the following SCE for this component: 



+? C(2) ^ E tl W ^ - ki) W^ 2 (ki), (6.20) 

where the parameters and depend on A and G\ (see App. [S]). On the basis of the transformation properties 
of this equation under common rotations of the replicas (and bearing in mind the coupling between positional and 
orientational degrees of freedom) , we propose that the solution has the form 

W^ZzM = (-l) mi « mi+ma ,ot« (1) (k 2 )+C 1 **r ^ (2) (k 2 ), (6-21) 

parametrized in terms of the two as-yet unknown functions w^- 1 ' and w^ 2 \ which each depend only on k 2 . By inserting 
this proposed form into Eq. ( 5.20| ), contracting (on the indices m\ and 7712), first with (— l) mi <5 mi + m2 ,o an d then with 
fc™* fc^ 2 2 * , and considering the limit k Q = k Q — > (with k 2 fixed and arbitrary), and retaining only terms that 

contribute to the leading-order value of Wi^iim (k), we arrive at the pair of coupled (Type II inhomogeneous) 
integral equations 
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c( 4 V 2 >(k 2 H-i^(k), 



c (4) - (1) (k 2 ) = ^ E kl ^( k - fci) ™ (1) (k?) + ? ^{ E £l " ki) W^,(ki)} 



+9 



E tl fci a ^A^> (2) (kDW(k-ki) 



(6.22a) 



(6.22b) 



where the parameter C^ 4 ) depends on A and C\ (see App. |5|). The symbol {• • -} g denotes the result of extracting 
from the quantity inside the braces the coefficient of the term proportional to the isotropic rank-2 spherical tensor 
(— f )" ll <5 mi +m 2 ,o- [To extract this part, we take the limit k" — k a ^ with k 2 fixed and arbitrary, and contract 

\mi. 



with (— l) mi <5„ ll+TO2 .o/3.] To find the leading contributions to and we first read off the value of the latter 
from Eq. ( 3.22a| ). We then use this result to eliminate from Eq. (6.22t), observing that we can omit the first 
term on the right hand side of Eq. (6.22h), and perform the remaining summations to arrive at the results 



w ( 2 )(k 2 ) 
™«(k 2 ) 

9k{6) 



3 v^F 



C (6) £ 2 



<Wir(0) e" k /e9 , 

_ , d0/c(0)e- fc2/ee , where 
47T 4 J 

de 1 7r(9 1 )de 2 TT{9 2 )e 1 9 2 s(9-(e 1 



i)) = ((M*))*(M*)))(*). 



(6.23a) 
(6.23b) 
(6.23c) 



where the symbol * represents Laplace convolution. Note that, as can be anticipated from our perturbative approach 
to obtaining the anisotropic components, the results for and are constructed from the universal scaling 
function 7r(0). 



3. Self- consistency equation: General case 



Having obtained the isotropic component of W, as well as the two lowest angular momentum anisotropic com- 
ponents, we now address the task of establishing the general form of W, along with an algorithm for obtaining the 
leading-order contributions to W for arbitrary values of its arguments. We begin by proposing the following structure 
for the general form of the leading-order contribution to W(k;£,m): 



W(k;£,m) 



Here Vs A (k) is a certain homogeneous polynomial in km in which all terms are of order £ in k. (We remind the 

reader that km stands for y/ (2£ + \)/A-k k l Yi m (k/fc), and is thus of order £ in fc.)In addition, 7^ m is a certain 

rotationally-invariant tensor depending only on £ and rh, and i^i m {9) is a distribution. All the unknown ingredients 

[Pl m (k), T i jh , and 7^(0)] will be obtained below. 

To illustrate these notions with a concrete example we note that W"^ " m (k), constructed in Sec. VIE 2 has such 
a form, with 



d9-K lm {9) 



(6.24) 



^>(k)ll 



£=(0,. .. ,0,1,1,0,. ..,0) 

f=(o,.. .,o,i,i,o,. ..,o) 

7r /,Th(^)li=(0 1 ...,0 1 l,l,0,...,0) 



lr(5) k a 



3 

H 



1 * I.OL2 * 
mi "'7712 ' 



1 A 



k{6), 



4V47T 



(6.25a) 
(6.25b) 

(6.25c) 
Note that 



where £ = (0, . . . , 0, 1, 1 , 0, . . . , 0) indicates that £ ai = £° 2 = 1, with £ a = in all other replicas 
Vj; ^j(k)||_( oiio o) ^ s °^ or< ier £ = 2 in components of k. 

To show that the leading-order contribution to the general component W{k;£,m) in the vicinity of the amorphous 
solidification transition does indeed have the form ( 6.24 ), for all values oi£, rh, and k, we proceed by full mathematical 
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induction on the multi-index £ jn|. We note that VF(k) has this form with ■ppg(k) = 1 and 7TQg(0) = 0, thus 



establishing the base of induction. To establish the step of induction we assume that W(k;£,m) has the form (6.24) 
for all values £ < £o (by which we mean £ a < £q for all a, and £ a < £q for at least one a). We then examine the SCE 
for W(k ;£n,m ): 

yn 

ki 

1 /4tt^ 2 



kiii,rfii oo,oi=l 



,2N 2 



2 V V" 



q J2 E W^Cki^^mOWCkaJ^,^) 

ki ,£i ,rhi k 2 ,f2,«i2 

A 

1 

X 



-1 ^ ( e -*0.6^s o Je^^ (626) 



ao,ai,a2 — 1 



As we are are only concerned with leading-order contributions to W, we truncate the £ sums in Eq. (3.26) so as to 



include only terms with I < Iq. We now examine the contributions to W^(ko, £q, mo) coming from linear couplings to 



lower-angular momentum components of W [i.e., the terms in the right-hand side of Eq. (6.26) that are linear in W 



Consider the coupling to W(k,£,rh). By the translational invariance of the correlator, performing the c average in 



Eq. (6.26) we establish that k = ko. The remaining correlator factorizcs on the replica index, becoming 

n 

II (IW-S.) W-s^e^) (fl.27) 

By introducing the Rayleigh-expansion, Eq. (|Al]) , for the factor exp ik$ • s™ we see that this correlator is nonzero 
only for the terms of order £ la in the Rayleigh expansion for which the angular momenta £' a , £q, and £ a can sum to 
angular momentum 0. As we are only interested in the leading-order contributions to W, and the £' a term in the 
Rayleigh expansion is of order (fe^) , we need only keep the lowest angular-momentum term (i.e., £' a — £q — £ a )- 
Thus, from each replica we shall pick up the factors (k§ )m^- m a, provided, of course, that |mg — m a \ < £q — £ a . 

Together from all replicas, these factors will give a factor, depending on the {km}, that will be of order Iq — £ in k. 
To leading order, the term that corresponds to the coupling of W(k ,£o,rho) to W(k,£,rh) will thus become, 



II (W%-^ C Z, mo ;l, m ) ^>( fc o) W(k ), (6.28) 

\a=0 



where the Cf - are constants, for which we cannot in general provide a closed-form result, but which we can, 
however, compute, should we decide to construct some component of the order parameter explicitly, as we have, e.g., 
done for ^1^; 

-A- s EIq (^o )m°-m° ^ s °^ or dcr £ — £ in ko, and as A (k ) is, according to the inductive assumption, of order £, we 

see that this contribution is of order £q in ko. Taking the term corresponding to £ — to, rh = rho over to the left-hand 
side, adding all remaining terms, and making the definition 

*W*o) = E E ( ft C U ^ (ko)^-2) r l7h (k ), (6.29) 

h \a=0 / 



we see that, to leading order in fco, the linear contribution to W(ko; £o, rho) is indeed of the form (6.24), with Vi Q ?fln (ko) 
being of the correct order in fco- As, due to the fact that q is small near the transition, the linear contribution is 
of lower order in e than the quadratic terms for values of k 2 of order e, we have thus also established a recursive 
algorithm for determining jha (ko). Illustrations of its use can be found in Apps. ^Jand^. 



We now examine the quadratic contribution to W(ko;£o,mo) in Eq. ( 3.24 ). This term is only the dominant one 
for k§ < e and, hence, we need only obtain it in the limit k^ — > 0, which is equivalent to extracting from it the 
leading-order contribution to the part proportional to the (k - independent) rotationally-invariant tensor 7| m . 

As we are only interested in the leading-order behavior of W, in the summations over £ in the quadratic term 
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of Eq. (6.26) we need only include the linear contributions to W(k;£, m) [i.e., for the purposes of computing the 
quadratic term of W(ko;£o,m ) we can set W(k;£, rh) = T>i A (k) W(k) for £ < £ ]. We now study the quadratic 
coupling of W-^ko; £o, rho) to lower angular-momentum components. In each replica a we have four sources of angular 
momentum: two {£" and 1%) coming from the two components of W [i.e., W(ki; £\, mi) and T / K(k2; £2, m^)]; and two 
(£'" and £'% ), one coming from each of the Rayleigh expansions of the "shift" factors exp iki • s ai and exp ik2 • s Q2 . As 
we are only interested in the leading-order behavior of W, we need only consider the case when £f +£% +£' 1 a +£ l 2 = £%■ 
Each source of angular momentum £ brings with it a factor of k e . Hence, multiplying all factors from all replicas 
together and assembling all terms in the summations we shall obtain 

^~"{ E fcl ,fc 2 FfriM) ^(ki) W(k 2 ) 5 fci+t2|to ^ 5 fc2>0 }^ (6.30) 

where the operation {• ■ -}r denotes the extraction, from • • •, of the part proportional to 7g Aq , and where T is a 

polynomial function of fci and k 2 in which all terms are of order £q. Provided, that £q and rho satisfy conditions for 
macroscopic rotational invariance (note, in particular, that MM requires £0 to be even and, hence, that the quadratic 
contribution vanishes for £ od d), we can perform the summations, indeed obtaining the claimed structure for the 
second-order term in Eq. ( |3.24| ). Thus, we have established a recursive procedure for obtaining tt^ A (#)• Keeping 

in mind that typical values of kp are of order e, that q — 2/3e, and that W(k) is of order unity we verify the scaling 
of the quadratic term to be e^°' 2 ' +1 , thus completing the step of induction, viz. the second and final element of our 
proof. 

In conclusion of our discussion of the solution of the order-parameter SCE we note, in passing, that all the assump- 



tions that we have made in Sec. VIA regarding the scaling (with e) of various quantities are verified, a posteriori, by 



the solution that we have obtained, thus establishing the self-consistency of these assumptions. 

VII. PHYSICAL INFORMATION ENCODED IN THE ORDER PARAMETER 

A. Introduction 

We have constructed the so lutio n of the order-parameter self-consistency equation in the vicinity of the amorphous 
solidification transition, Eq. ( |6.24 ), and have obtained explicit solutions for the two lowest angular- momentum com- 



ponents of the order parameter, Eqs. ( |6 .18 6.23a , 6~23~fc| ). We are thus in possession of a range of statistical inf ormat ion 



about the amorphous solid state, this information being encoded in the order parameter, as discussed in Sec. Ill A. In 
this section we will extract some of this information explicitly, from the two lowest angular-momentum components 
of the order parameter, provide a strategy for obtaining other statistical information about the system from higher 
angular-momentum components, and discuss the scaling of the order parameter with e near the transition, as well as 
the implications of this scaling. 

Our statistical diagnosis of the structure of the amorphous solid state in the vicinity of the amorphous solidification 
transition is made in terms of the moments of the joint distribution function P that collects together the localization 
characteristics of all the particles in the sample and their orbitals, averaged over realizations of the disorder. As we 
shall see, we are unable to construct the entire distribution function, or even to construct a closed-form expression 
for arbitrary moments. However, most of the useful information about the system can be obtained from low moments 
of the distribution, resulting from the two lowest angular-momentum components of the order parameter, which we 
have obtained explicitly. As for the information encoded in the higher angular-momentum components of the order 
parameter, we do not extract it explicitly. We do, however, describe the kind of information that could be obtained 
from them, as well as provide the general procedure for doing so. 

B. Encoded information: Isotropic sector 

We begin by discussing the statistical information, encoded in the isotropic component of the order parameter 



VF(k 2 ). In Sec. VIP we have obtained this component in the vicinity of the amorphous solidification transition and, 
consequently, the reduced distribution of (inverse square) localization lengths p(r) associated with localized particles. 
The fact that p(r) takes on a scaling form, and that the scaling function tt(9) has a well defined peak with location 
and width both of order unity, allows us to establish that (as e — > 0) r scales as e and, accordingly, £ scales as e -1 / 2 . 
For example, any reasonable choice for a characteristic value of £, say J °° drp(r) r -1 / 2 , scales as e" 1 / 2 . 
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In the present context, while the emergent distribution p(r) is found to have striking scaling properties, there is the 
suspicion that changes in the details of the model will lead to changes at least in the details of the scaling function ir(9) 
or, perhaps, in the scaling property itself. However, it has been found in the context of vulcanized macromolecular 
matter that the scaling property, as well as the precise form of the scaling function, are robust, universal features of 
the mean-field theory, verified by independent computer simulations. Moreover, this universality has been found to 
have its origins in the symmetries of the appropriate Landau free energy and the divergence (at the transition) of the 
characteristic localization length. 

Being in possession of the entire distribution of localization lengths provides us with a surprisingly rich character- 
ization of the positional aspects of the amorphous solid state. It is striking that the distribution is universal, not 
only across the macromolecular systems where it was first found, but also in the present setting of vitreous media. 
Although we have obtained the distribution via analysis of a specific semi-microscopic model of vitreous media, we 
anticipate that, here too, the result will have a broader domain of applicability. Moreover, given the emerging picture 
of orientational order as order slaved to the underlying positional order, it is not surprising — and indeed we shall 
demonstrate this point below — that all other statistical descriptors of the amorphous solid state are also constructed 
from the universal function tt(9). 



C. Encoded information: Low angular-momentum anisotropic sectors 



1. Angular localization 



We now discuss some of the specific physical information that can be obtained by examining our explici t solut ions 
for the two lowest angular-momentum components of the order parameter, which we have obtained in Sees. VIE 1 and 
VIE 2. The first piece of information concerns the an gular localization of the orbitals, without regard to the positional 
localization of the particles. As we recall from Sec. Ill A, such information is accessed via the order parameter for 
k = 0, and is described by the distribution of the collection of characteristics {XV m }. The lowest angular-momentum 



order-parameter component that provides access to this information is W" 



mi I1712 



(k), the solution for which is given 



by Eq. (5.21). Evaluating at k = gives 



QWZ 1 T m2 (0) = (-l) mi S mi+m2 ,o 



e 3 

4~7T 6 



d9n{9). 



(7.1) 



Now, recalling the interpretation of the order parameter developed in Sec. [II A, we arrive 



f 1 N 1 A 

q / drd{E}P{r}P(r,{I]},{r}) E ljmi E 1>ma = ^EtE^M) (C 2 fe« 

** A 1 1 



j=l o=l 



C (6) £ 3 poo 

(_l) m ^ mi+mai0 -=- / dS K (9). 



47T 6 J 



(7.2) 



Notice that this characteristic of the angular localization of the orbitals is essentially the order parameter traditionally 
used to describe the directional localization of magnetic moments in the spin-glass state. In fact, if we recall the 
spherical-tensor decomposition of the scalar product of two unit vectors, si -s 2 = ^ m (— l) m Yi m (si)Yi_ m (s 2 ), then by 
appropriately contracting Eq. ([7.2]) over mi and vn^ we obtain the familiar characterization of directional localization: 



1 N 1 A 



N ^ A- 

j=l 0=1 



3C< 6 ) e 3 

47T 6 



o19k{9). 



(7.3) 



2. Angle-position fluctuation correlations 



Further specific physical information concerns the degree to which the thermal fluctuations in the orbital orientations 
are correlated with the thermal fluctuations of the particle positions. As discussed in Sec. Ill A, to extract this 
information, which is encoded in the distribution of the collection of functions {r£ m (k)}, we examine the order- 
parameter component W" m (k) : 
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qW^ ni (k)=q J dTdmV{T}P(T,m,{T})T hmi (k^)e- k2 / 2 
1 N 1 

T- 

^ A 

3=1 a=l 

d6»7r(6»)e- k / e 



3=1 

2e iC« 



^- jk <n.(c 3 -^-) _ ^-^-.(c,-^))) (^(s^) - (y x * mi Ka)) 



U a l * 
"ill! 



'k,0 



3 \/3 V4vr 

By considering the derivative with respect to k Ql , taking the limit k > 0, and contracting, we obtain 



(7.4) 



1 N 1 A 

j=l a=l 



C (l) 



(7.5) 



characterizing the anticorrelation of the thermal orientation-fluctuations of orbitals with the thermal position- 
fluctuations of the particles to which the orbitals are attached. 



D. Encoded information: Higher angular- momentum anisotropic sectors 



In the previous two subsections we have obtained explicit physical information about the amorphous solid state 
in the vicinity of the solidification transition from the three lowest angular-momentum components of the order 
parameter. We now discuss the generalization of the ideas that we have used in extracting that information. 

Take, for instance, the task of obtaining information about the orientational localization of the orbitals. Although 
the spin glass order-parameter-like quantity [(s) • (s)], constructed from [(Yi mi (s)) (yi m2 (s))] , provides the simplest 
characterization of the orientational localization of the orbitals, one could also study its higher angular momentum 
analogues, such as, e.g., [(l^mi (s)) (l2m 2 ( s ))] or, in general, [Yl a (Ye<x m <*(s)) ] . Such quantities could be useful if, e.g., 
one wished to study the properties of the formed covalent bonds separately from those of the unbonded orbitals. Such 
quantities are obtained as appropriate moments of the joint probability distribution P which, in turn, are extracted 
fro m the c orresponding components of the order parameter. Specifically, W(0;£, rh) yields [J7 (Y£a ma (s)) ] . In 
Sec. VI E 3 we have established an algorithm that allows us to obtain the leading-order contributions to all components 
of the order parameter. Accordingly, we also have an algorithm for obtaining a large class jl6) of the moments of the 
joint probability distribution P, including [J\ a (Y^a m c.(s))] for arbitrary I and rh, via 



j—l a—l a—0 



= q / dTd{X}V{T}P(T,{i:} > {T})£ i .=qW{d;i 



(7.6) 



where was used to denote Y[ a ^i a m a - 

Similarly, higher-order characterizations of the thermal orientation-position correlations, i.e., the moments 



q \ drd{E}I?{r}P(r,{S},{r})f i .(k)exp(-k 2 /2r) 



(7.7) 



can be obtained from q W(k; £, rh). [Note that f/~ (k) denotes Y[ a ^e a m a (k a ).] For £ odd these moments are simply 
equal to q W(k; £, rh). For £ even, they are equal to q W(k; £, rh) — q W(0; £, rh). Then, to find such characterizations, 
just as we did in the previous section, one needs to extract the coefficient of (k a Y ma , and pass to the limit k — > 6. 
For example, by applying this procedure to lm s ) we wou ld obtain 



1 N 1 A 



■in 



(7.; 
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E. Scaling and its indications 



We now address the manner in which various physical quantities reflecting the orientational ord ering sc ale with e. 
To do this we examine the scaling of various components of the order parameter, established in Sec. VI E 3 . Hence we 
have the following scaling for various moments |L6| of the joint probability distribution P associated with localized 
particles (£ ^ 0): 



dr d{£} V{T} P(t, {£}, {T}) £*. ~ e 



dr d{£} V{T] P(t, {£}, {T}) f (k) exp -k 2 /2r 



l+t/2 



(7.9a) 
(7.9b) 



the latter being valid for k ~ e. Note the exponent 1 + 1/2 in Eq. (7.9a) 



We now describe a plausible physical scenario that is consistent with this scaling of P. Near to the transition 
the network has many long chains consisting of twice-bonded particles, with only the occasional more highly-bonded 
particles linking them. Most of the localized particles are on extremely mobile segments. Orbitals attached to particles 
on such segments are likely to be even less localized orientationally than orbitals on less mobile segments, such as 
those near junctions between chains. Consider, for the sake of illustration, a dangling chain (i.e., one attached to 
the network only at one end). The orientational fluctuations of successive orbitals compound the fluctuations of 
orbitals further away from the junction, causing them to be successively more delocalized. This compounding of 
fluctuations so heavily suppresses orientational localization that it causes the scaling (and not just the numerical 
value) of the orientational localization characteristics E^ m to vary according to location in the network. Moreover, 
because orientational localization is heavily suppressed for such a large fraction of orbitals, the leading-order scaling 
of the moments of E^ m is dominated by contributions from the small fraction of better-localized orbitals (e.g., those 
near chain junctions), and is blind to the much larger fraction of less well-localized orbitals (e.g., those far from chain 
junctions). This partitioning of localized orbitals into better and less well fractions yields a picture consistent with our 
results provided we assume that the better localized variety co nstitu te a fraction of order e of the localized orbitals. 
This fraction manifests itself in the exponent 1+1/2 in Eq. ( 7.9a ) as the additional 1. This picture allows us to 
identify the following appealing scaling for the localization characteristics of the better localized orbitals: 



r, m (k) ~ ^ 



(7.10a) 
(7.10b) 



the latter being valid for k 2 ~ e. Whilst we cannot be certain that this scenario is a necessary consequence of our 
results, it is both consistent with them and physically plausible. 



VIII. CONCLUDING REMARKS 



The primary result of this Paper is a statistical characterization of the structure and heterogeneity of the equilibrium 
amorphous solid state that emerges due to the random permanent covalent bonding of the constituent particles. This 
statistical characterization takes the form of a joint probability distribution that ascertains the likelihood of finding 
a particle: (i) to be localized, (ii) to have a certain positional localization length, (iii) for a bond connected to this 
particle to have a certain orientational localization characteristics, and (iv) for the correlations between the thermal 
fluctuations in particle position and orbital orientation to have certain characteristics. 

We expect the emerging statistical characterization to be valid beyond the context of the model used to determine 
it. The reason for this is that, apart from certain simple dependences on physical parameters describing the particles 
in the network, this characterization is in fact a consequence of the order parameter that we have considered, the 
symmetries of the Landau-type free energy associated with this order parameter, and a limited number of further 
assumptions. 

We have constructed an analytical approach to the equilibrium structure of the amorphous solid state of a class of 
materials, such as silica gels, formed by the permanent random covalent bonding of atoms or small molecules. However, 
the accuracy and scope of our results is limited in the following significant ways. We have focused on equilibrium 
structural properties, we have worked near to the amorphous solidification transition, and we have computed within the 
framework of a mean-field approximation. It would be interesting to have a better understanding of the implications 
of these limitations, and to be able to obtain results beyond them. 
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APPENDIX A: SOME USEFUL MATHEMATICAL INGREDIENTS 



Owing to the presence of the factor exp(ik • s) in the effective Hamiltonian, we shall need the Rayleigh expansion 
for plane waves: 

oo e 

exp (tk • r) = 4tt^ £ i e j £ (kr) Y? m (k/fc) Y tm (r/r), (Al) 
e=o m=-i 



where jt{p) is a spherical Bessel function, i.e., je(p) = y/ir/2pJ l \ (z). The small- z asymptotics of jt(z) for small 



't+i 

are: 



j (z)~l-z 2 /6, (A2a) 
h(z)~z/3. (A2b) 



We also make use of angular momentum one spherical components of vectors: 

k* m ^^kY* m (Vk). (A3) 

Lastly, we need some special values of Wigner 3 — j symbols: if l\ + £2 + £3 is odd, then 

«M=0, (> 1 °) =i „ 1+ „„ o( - ir 1. (A4) 
/ \mi m 2 J V3 



APPENDIX B: AVERAGING OVER THE DEAM-EDWARDS DISTRIBUTION 

In this Appendix, following the discussion in Ref . JhJ , we perform the disorder-averaging of the n-fold replicated 
partition function Z n with respect to the Deam-Edwards distribution Vm{C) characterizing the distribution of the 



quenched random constraints C = {j e , j' e ', a e , a' e j M _ 1 (i.e., bonds), defined in Sec. [V B . We start from the definition 



of the disorder average as a weighted average over all possible realizations of the disorder: 

[Z n (C)]=Ti- c V M (C)Z n (C), (Bl) 
where Trc denotes a trace over all possible realizations of the disorder C, defined via: 



M 



^ -4EnlEE )•••• (B2) 



A/=0 e= 



Here, A^ is a normalization factor, which we shall compute below, and the sum over M collects contributions from all 
possible numbers of permanent covalent bonds. Note that our definition of the Deam-Edwards distribution Vm{C) 
includes the weig hting factor (2irp 2 V/N) M /M\ which is a quasi-Poissonian probability distribution for the number of 
bonds M, controlled by the parameter p 2 . Hence, we have 

M=0 e=l \ j e ,a a jl,ai / e=l 

(B3) 
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It is worth noting that the sums over j e ,j' e , a e , and a' e are not restricted to preclude the situation where a particular 
orbital a on a particular particle j is bonded to two (or more) other orbitals attached to other particles, an unphysical 
situation, as covalent bonds are typically regarded as "saturable" i.e., any orbital can be bonded to at most one other 
orbital. However, near the amorphous solidification transition the bond density is low, thus we expect the occurrence 
of such multiply-bonded orbitals to be negligibly rare and thus not having much impact on the theory. 

Introducing replicas, denoting by ( • • • } N 1 the thermal average weighted by a Hamiltonian that does not couple 



the replicas, and resumming the power series in M, we can rewrite Eq. (B3) in the final form 



[Z n (C)] =Z n+1 /Z x , where (B4) 



55 ( exp (w £ £ II 5(3) K + K. - c r ~ K>0 A(2) (s" Q - ) ) ■ (B6) 

\ \ i,i'=la,»'=la=0 / / N,n+1 



N 



APPENDIX C: DENOMINATOR FOR THE SELF-CONSISTENT EQUATION 

In this Appendix we evaluate the denominator of the right-hand side of the SCE for n(k;l, m), Eq. ( |6.2| ), namely, 
the quantity 

( ex P E ^(k; *> rh) i 2 «'' k s > „„■: §«)) \ . (CI) 

To do this, we expand the exponential in a power series and consider the r-th order term: 

l^ (^H 2 ) r q r V- nr E ■■■ 2 ^(ki;^,m 1 )I^(k 2 ;£ 2 >2) ■■■VK(k r ;i r ,m r )e lil - e e^ 2 - 6 ...e 4 ^- 6 

* kiiimi 142^2^2 k r i? r m r 

X ^ E eitl - S01e * 2 ' Sa2 --- eifcr ^^ 1 (- § aJ---% Af .(-§a r )) • (C2) 

Ol,...,Or / 1)7l+ i 



Noticing that the quantity (|C2j) factorizes over the replicas, and passing to the replica limit n — > 0, we rewrite 
(i(47r M 2 ) r g r ^ E ■■• E ^(k 1 ;£ 1 ,m 1 )^(k 2 ;^ 2 ,m 2 ) ...iy(k r ;€ r ,m r )e 4kl - c e 4k2 - c ...e ik - c 



ki^mi k2^2^2 k r / r m r 



x E e ikl - s ^e ik2 - s ^...e ik -^^ imi (-s 01 )...^ rmr (-s ar )\ . (C3) 

ai,...,o,. / 1 j 

By the MTI of T4^(k; m) we know that k^ = 0, which means that k^ = 0. Also, using MRI of W(k; £, m) we see that 
if k = 6 then ( |C3| ) can only be nonzero if I = 6 also. Knowing that W(6; 6, 6) = 1/V47T, we get, for the r-th order 
contribution (/x 2 ) q r /r\. Resumming the power series we obtain, finally 

( ex p E w & ™) \ E elk '^ («.)) ) = E ^~~r~ - -p(/^)- 



APPENDIX D: DEFINITIONS OF THE C (JV) CONSTANTS 

The constants are defined as follows: 
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C (0) 


= A-^l + ^-lJCi) 


(Dla) 


C (l) 


1 

~ i + c(°) 




C (2) 


= A. fi + (A - l)Ci 2 ) 


(Die) 


C (3) 


= C (2) _ 2C (2) C (1) 


(Did) 


C (4) 


1 - 
~ CP) 


(Die) 


C (5) 


C w 


(Dlf) 


C (6) 




(Dig) 


~ 12C( 2 )CW 2 ' 8CW 



They encode information about the strength of the mutual repulsion of orbitals. 



APPENDIX E: ANISOTROPIC SELF-CONSISTENT EQUATION: ANGULAR MOMENTUM ONE IN 

ONE REPLICA CHANNEL 

In this Appendix we study in detail the SCE for the W"^ (k) component. We start with the full form: 



x 4 E (e-^nWC))^^)^ 1 -^* 1 ^^^) 



ao,ai — 1 \ q^qi 



l,n+l 



+ K^") ^ ? ? ^(ki;ii>i)^(k 2 ;i 2 > 2 ) 

kl ,£l ,7TM k 2 ,<?2 ,™2 

E ( e "* £ ft (^oo(s-)) ^(CK^e* 1 - 



Oo,Oi,a2 = l \ Q^Ql 

x%(§a 1 )e ik - e e^'^%(s 02 )\ (El) 

/ l,n+l 

We now proceed to study this equation order by order, starting with the linear term. As we are not interested in 
feedback on this component from higher angular momentum components we may truncate the angular momentum 
sums to include only terms of £ 2 < 2. Thus we must include couplings to the components W"^ (k) and W(k) only. 
We study the coupling to T4 / (k) first (i.e., l\ = rh\ = 0). Clearly then in replica oei we must expand the shift-factor 
cxp(«k" • s^) to angular momentum 1 and put this factor equal to unity in all other replicas, obtaining 

ki 

1 A I . . / 1 \ n / 1 \ 1 /k c 



E (e-*^^=j Y? mM) [ ) E ^(s-)^, 

o ,oi = l \ m'=-l x 1 7 / i n +i 

= (l-g) M 2 z^^(k)^ ^ (^ + (1-^)^)^^^) 

ao,ai = l 

= -j=C(°)Ty(k)fc« 1 1 *- (E2) 
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We now study the linear term corresponding to the coupling of W"^ (k) to (k). When a\ ^ a-i this term, 

being of order k 2 , is subdominant. Thus we need only consider the term in the angular momentum sum where £ ai = 1, 
l\ being in a ll oth er replicas. We examine this term in the sum, making use of the value of the two-orbital correlator, 
given by Eq. (6.17): 



ao,ai — 1 
A 



-ik-c ^iki-c 



Y* (<i ai \ 
1 lmi V B a I 



1 



4tt/ / ln+ i 



ffQ ,0-1 — 1 



= -(i- 9 )c(°) A1 2 vy- i (k). 



(E3) 



We next proceed to examine the second-order term in Eq. (El). Again, we truncate the angular momentum 
sums. As we are only interested in the leading-order contributions to W"^ (k) there are only two cases to consider: 



l\ = £2 = 6 and l\ = 2, l\ = 0. The former case, however, results in a term that is subdominant to term (E2). We 
thus consider the latter case: 



1 



/4V X 2 



W{^)W^(%) 



ki,k 2 



A 2 E 



ao,a± — 1 



1 / 1 



^lmi( S cto) ( n — ) s ai) 



1,71+1 



(l-^^E^^^^k-kx)-! £ (*«,« 1 + (l-*-,.x)^i 



ao ,ai — 1 



(1 - g)C<°> ^ Yl ^( k i) W ^ ( fc - ki). 



(E4) 



Keeping in mind that due to the symmetry l\ <-> I2 we must double the term (E4), we now assemble all the terms, 
dropping sub-leading contributions: 



qW«\ n (k) = C<°) ( - gW^ (k) + ^=gW&) fc** ~ <? 2 ^ E kl - kl ) W ^ ( kl )) • ( E5 ) 
Cancelling a factor of q leaves us with Eq. ( |6. 16 ) . 



APPENDIX F: ANISOTROPIC SELF-CONSISTENT EQUATION: ANGULAR MOMENTUM ONE IN 

TWO REPLICA CHANNELS 

In this Appendix we study in detail the SCE for the W / im 1 im 2 (k) component. We start with the full form: 



- £ (e^' 6 II (^o(s«J)r; mi (s-)r; m2 (s-)e^- 6 e^-- 1 f irii (s ai ) 



1 / V N ' 



l,n+l 



2(l^) q2 E E ^CMi^i)^;^) 



ki,^i,?fii k 2 ,£ 2 ,m 2 



A 3 



ao ,ai ,02 — 1 
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ik 2 c ik 2 s a 



(Fl) 



l,n+l 



As we did before in App. ||] for the SCE for * (k) we proceed to study the linear terms of this equation. Again 
we truncate to include only the components of W of angular momentum smaller than the angular momentum of 
Wimiim 2 (k)» i.e. £ a = for a f a u a 2 and£ Q \£ Q2 < 1. Thus we find couplings to terms W(k), (k), W^k), 

and Wimi im 2 (^-) • Expanding the shift-factor exp(ik" -s^) to angular momentum or 1 as appropriate in each 



repli ca, pro ceeding as in Eqs. (E3) and QE2Q, and, finally, inserting the value for W^^k) that we obtained in 
Sec. VIE 1, we obtain the linear contribution (terms are in the order listed above): 




— tk-c iki-c 



^lmi ( S ao )^lm 2 ( S oo ) X 



(4™ f j t (k^ ) ji (k^ ) £ F M (C ) l? m . 



ki c 



k* 1 



Ylm' 2 ( s af ) ^lr 



+^ 1 (ki) 



1 



+^ 2 (ki) 

+^i^(ki) 
1 



4tt 
1 



l 

Y lmi (-s%)4ni ji (k^) Y, y i-'(O y i*, 

m'=— 1 



Y lro2 (- S S)4 7 rzj 1 (fc 1 ai ) E ^'K 1 )^; 

m'= — 1 



ki c 



k^ 



kl c 



A-" 1 



ilmi( s oi)^lm 2 ( s ai) 



l,n+l 



(F2) 



We next proceed to study the quadratic contribution to (k) from Eq. (Fl). Truncating the angular 

momentum sums as usual we see that there are the following four cases to consider: 

(i) l\ = £2 = 0, corresponding to the coupling to W(kx)W(k2) 

(ii) £\ =0,11 = 2, corresponding to the coupling to W(ki)W"^ li (k 2 ) 



(iii) £\ 



P 2 

l 2 



2, corresponding to the coupling to PF" I 1 „ 1 (ki)W / ™^ l2 (k 2 ) 



(iv) £1 = QJl = 4 corresponding to the coupling to W(ki ) W?* „ (k 2 ) 
We study these terms in more detail, as in Appendix ||| Eq. (E4). Term (i) is clearly subdominant compared with 
the second term in Eq. (F2). Writing out the remaining terms: 



term(ii) 



x/3 2 



(1 - q )c^ iL E ? 2 C 2 *^(ki) w? mx (k - k x ), 



(F3) 



term(iii) -.1^(1- q ) C ^^ 



2 

4tt 



V" ^ki 

,V(2) M 



term(iv) :^(1 - g^C^ ^ ^ 9 2 W(kl) W&^k - k x ). 



(F4) 
(F5) 



Term (ii) does not contribute to the leading-order behavior of W"^ ni (k) as for k 2 of order e or greater all quadratic 
terms are subleading compared to linear terms, and for k 2 <C e it will, because of the k^ 2 * factor, be subdominant 
compared to terms (iii) and (iv) . Keeping in mind that because of the symmetry t\ <-> £2 we must double the term 
(iv), we reassemble the pieces of the SCE for W"^ 

llmst )' dropping subleading contributions: 



An 



yn Z-^i 

ki 
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(F6) 



ki 



Cancelling a factor of q leaves us with Eq. (6.20). 
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